Relativistic Mott criticality in graphene 



0^ ■ 

o : 

o , 
(N . 

D ■ 
00 ; 



I 

o 
o 



> 
0^ 



O 

On 
O 



Igor F. Herbut^, Vladimir Juricic^, and Oskar Vafek^ 
^ Department of Physics, Simon Eraser University, Burnaby, British Columbia, Canada V5A 1S6 
^ National High Magnetic Field Laboratory and Department of Physics, 
Florida State University, Tallahassee, Florida 32306, USA 

We formulate the effective Gross-Neveu- Yukawa theory of the semimetal-insulator transitions on 
the honeycomb lattice and compute its quantum critical behavior near three (spatial) dimensions. 
We find that at the critical point Dirac fermions do not survive as coherent excitations and that 
the ~ 1/r tail of the weak Coulomb interaction is an irrelevant coupling. The emergent Lorentz 
invariance near criticality implies a universal ratio of the low-temperature specific heats of the 
metallic and the rotational-symmetry-broken insulating phase. 



I. INTRODUCTION 

It is well known that the large overlap between the Pz 
orbitals from neighboring carbon atoms makes graphene 
an excellent conductor.^ Nevertheless, one can conceive of 
situations in which the relative strength of the repulsive 
Coulomb interaction between electrons would be higher 
so that graphene would turn into a Mott insulator. A 
recent calculation suggests that just taking the graphene 
sheet away from the substrate may gap out the Dirac 
pointSi-^? If so, varying the dielectric constant of the sur- 
rounding medium over a sufficiently wide range could 
in principle be used to tune through the metal-insulator 
(MI) transition in graphene. Another possibility would 
be stretching the sheet to reduce the hopping between 
the Pz orbitals. Such a quantum phase transition would 
be the analog of the Higgs mechanism for gauge-neutral 
fermions in particle physics, with the Higgs boson here 
as a composite field. It would represent maybe the sim- 
plest example of fcrmionic quantum criticality, in which 
the gapless fermions exist only near the isolated points 
in the momentum space. A finite gap would also make 
graphene more interesting for potential applications in 
electronics.* 

A unique feature of the Mott phase in graphene is that 
it may come in several varieties: from the familiar Neel 
and staggered-density phases^^ to the more exotic in- 
sulators that break the time-reversal symmetry (TRS)ji2i 
There have been several studies of the MI transition in 
graphene i^ii'^'^iii^'iiii^ Nevertheless, several fundamen- 
tal questions still await answers. Among these, the fol- 
lowing qualitative issues seem particularly pertinent: 1) 
what is the role of Dirac fermions in the critical behavior, 
2) how does the criticality depend on the nature of the 
order parameter (OP) in the Mott insulator, 3) what is 
the fate of fermions near the critical point, 4) is the long- 
range tail of the Coulomb interaction relevant? Ques- 
tions 1) and 3) echo some of the central themes of the 
wider field of quantum critical phenomena j^'^i^* whereas 
the last question, as we will argue, is related to the clas- 
sic problem of triviality of the continuum limit in non- 
asymptotically free field theories. 

In this paper we present an effective theory of the 
Mott transition in graphene which allows us to address 



these and related issues in a controlled and transpar- 
ent way. Our action contains both the self-interacting 
bosonic (or "Higgs") OPs and the Dirac fermions, cou- 
pled by Yukawa-like terms. It represents a simple modi- 
fication of the Gross-Neveu theory, derived previously in 
the large-N limit, ^ but with a crucial novel feature: there 
is an upper critical (space) dimension in the problem of 
three. This allows one to perform the e = 3 — d expansion, 
with the Higgs and the fermionic fields at all stages of the 
calculation treated on the same footing, thus placing the 
Mott criticality in graphene at the same level of rigor as 
the textbook 'I'"' theory. We find the MI transition in 
graphene to be of the second order, and to be governed 
by the critical point laying at a finite Yukawa coupling 
(Fig. 1). Although crucial for the critical behavior, Dirac 
fermions acquire a small positive anomalous dimension, 
so that the residue of the quasiparticle pole continuously 
vanishes as the transition is approached from the metallic 
side. Whereas the transition may be tuned by increas- 
ing the strength of Coulomb repulsion, its ~ 1/r tail is 
in fact an irrelevant perturbation to the leading order 
in e. We determine the dependence of the critical expo- 
nents on the broken symmetry of the Mott phase. Our 
analytical results compare favorably with those of a re- 
cent numerical worki^ The emergent Lorentz symmetry 
near criticality implies the existence of a new universal 
quantity: the ratio of the low-temperature specific heats 
of the semimetal and of the rotationally non-invariant 
insulator. Finally, possible analogies between the Mott 
criticality in graphene and the chiral symmetry breaking 
in 3 -f- 1-dimensional quantum electrodynamics {QED^) 
are noted. 

The remainder of the paper is organized as follows. In 
the next section we discuss the order parameters and the 
effective field theory near the transition. In sec. HI we 
perform the analysis of the problem near 3-f 1 dimensions, 
and in sec. IV we discuss our results, as well as the 
universal amplitudes specific to the problem at hand. 



II. EFFECTIVE THEORY 

Let us begin by fixing the necessary notation. Intro- 
duce the Dirac fermion for spin- 1/2 electrons as = 



2 



(^'1,^1), where each spin component is defined as 




FIG. 1: Schematic flow in the critical plane for Mott transi- 
tions in graphene. For any positive bosonic quartic coupling 
A the transition is continuous and governed by the fermionic 
critical point F. g is the Yukawa coupling between the OP 
and Dirac fermions. The other marked fixed points are the 
Gaussian (G) and the Wilson-Fisher (WF). The (unmarked) 
bicritical fixed point is unphysical, as argued below. 

where K = {1, 1/^/3) {2'K/a\^) is the Dirac point, and 
u and V are the Grassmann fields on the two sublattices 
of the honeycomb lattice. The reference frame is rotated 
so that qx = q ■ K / K. The non- interacting Lagrangian 
is then Lj — 4*7^ 9^4', with the 7-matrices defined as 
7o = -^2 <Xi (Tz, 71 = CTz (El (Ty, and 72 = ^2 <Xi cr^, and 
\]/ = v]/t^Q_ The ultraviolet cutoff is A ^ 1/a, where a 
is the lattice spacing. We have set the Fermi velocity to 
unity. 

In 2-1-1 dimensions there are two qualitatively dif- 
ferent ways to gap out the Dirac fermions.^ The first 



i'=<t>,X 



and the Yukawa terms that couple bosonic and fermionic 
fields are 

^y= {g^As^M^^^ + g^^t-^t ■ ^M^CT*}, (5) 



, vl{,K + q,ujn), ul{~K + q,uJn),vl{-K + g, w„)), (1) 

I 

one preserves the chiral symmetry (CS) generated by 
{73,75,735}, 735 = «7375, but breaks the TRS for each 
spin component separately. The time-reversal in the 
above representation is defined as ^'o- — > /t^'o-, where 
It = [px ® I2)K and K is the complex conjugation.^^ 
Choosing 73 = CTj. (g) ay and 75 = ay ® ay, the TRS- 
breaking OP may be written as 

with (ps and ipt as its (real) singlet and triplet compo- 
nents. The Pauli matrices act on spin, and 7-matrices 
act on Dirac indices. A lattice realization of has been 
offered long ago in terms of circulating currents between 
the neighboring sites on the same sublattice.J^ It was 
also shown recently that ^ is favored by the second 
nearest-neighbor repulsion on honeycomb latticeJ^ 

The second order parameter preserves the TRS, but 
breaks the CS: 

X=(x.,Xt) = ((**), (3) 

In the representation above, for example, the singlet cor- 
responds to staggered density, and the triplet to stag- 
gered magnetization. A finite Xs {xt) may then be in- 
duced by a large nearest-neighbor (on-site) repulsion;^ 
In this case a finite singlet component breaks the Ising 
symmetry of the sublattice-exchange, whereas the triplet 
breaks also the (spin) rotational symmetry. 

We may promote the above fermion bilinears into 
separate dynamical fields and write the effective ac- 
tion near the Mott transition as 5* = / drdxL, with 
L = Lf + Ly + Lh + Lc, where the bosonic part is 



tV^ -I- t^^t)i^t] + A^0,sV's + Av-,t (V^t ■ i^tf + A^,,tV'i??}, (4) 
I 

where the matrices = 1, and = 735. Note that L 
lacks the Lorentz symmetry due to the assumed generic 
bosonic velocities v^,s ^ w^,t 7^ 1- The form of L is dic- 
tated by the spin-rotational, time-reversal, and sublat- 
tice exchange (Ising) symmetry of the original problem 
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on honeycomb lattice. 

One may also account for the long-range tail of the 
Coulomb interaction by including the term 



(6) 



We allow for a general spatial dimension d, but in a man- 
ner that assures that the integration over the scalar gauge 
field a would introduce the ~ e^/r density-density inter- 
action between fermions in any dJ^ 

Integrating out the Higgs OPs in the symmetric phase 
t^,s,tip,t > would produce local quartic terms for the 
remaining fermionic degrees of freedom, and reduce the 
theory to the Gross-Neveu form derived before in the 
large-N limit The advantage of the Yukawa form is that 
it becomes renormalizable in 3-|-l dimensions, where both 
the Yukawa and the Higgs self-interaction couplings be- 
come dimensionless. The relationship between the Gross- 
Neveu and the Yukawa descriptions mirrors the one be- 
tween the non-linear and linear sigma-models in standard 
critical phenomena, and there is a substantial evidence 
that both describe the same critical point iS2, In the re- 
maining part of the paper we exploit this equivalence 
and study the effective Yukawa theory in d = 3 — e di- 
mensions, with an eye on e = 1 relevant to graphene. 

A possible objection to the pursuit of the above strat- 
egy may be that the Higgs bosons that involve 735 do 
not gap the Dirac fermions in 3 -I- 1 dimensions, where 
the CS is reduced to J7(l). If we were to compute the 
beta-functions directly in 2 -f 1 dimensions, on the other 
hand, at least perturbatively there can be no difference 
between those with x and cj) labels. This is because the 
matrix 735 commutes with all the 7-matrices which ap- 
pear in the fermion propagators.^ This leads us to con- 
jecture that the critical behaviors at x and (p transitions 
that share the same spin-symmetry are identical. 



III. e - EXPANSION 

With the above rational in mind, hereafter we consider 
only the CS breaking transition. Anticipating some of the 
results we set v^^s — %,t — 1 and = 0, and integrate 
out both the fermionic and the bosonic modes within the 
four-momentum shell A/b < (w^ -I- P)^/^ < A. There 
are two critical points, distinguished by the rotational 
symmetry of the OP: 1) at i^.^ = {t^^t = 00), and 2) at 
^x,t — i^x.s = 00)1^ They may be tuned by increasing 
the nearest-neighbor and the on-site repulsive terms in 
the Hubbard model, respectivelyj^ To one-loop order we 
find 



dgl 

dlnb 



dXu 
dlnb 



Xu{e-8gl)-Ai9 + S)Xl + 2gt 



(7) 



(8) 



where g„ = (?;^,„/(87r^A^) and A„ = A^,„/(87r^A'^), u = 
s,t, and S = {S = 2) for singlet (triplet) couphngs. 
The purely bosonic Wilson-Fisher fixed point is unstable 
in the g-direction, and the critical point lies at a finite 
Yukawa and quartic couplings (Fig. 1). The flows in both 
cases resemble the one of superconductors,^^ except that 
the bicritical point now lies in the unphysical region, so 
that both the singlet and the triplet transition are always 
second-order. 

Let us examine the stability of the critical point with 
respect to perturbations. First, breaking the Lorentz 
symmetry by a small difference in the fermionic and 
bosonic velocities is irrelevant: for 6u — I 
to the leading order in e we find 



« 1, 



4e 



dS„ 



d\nb S-7 



(9) 



Weak Coulomb interaction is also readily found to be 
(marginally) irrelevant 



de^ 
dlnb 



(10) 



similarly as in the purely bosonic theoriesJ^ Here we 
rescaled c^/Stt^ e^. The first term in the bracket de- 
rives from the usual polarization, and it is present only 
in d = 3, where the second term in Eq. (6) becomes an- 
alytic in momentumfi^ The remaining term comes from 
the renormalization of the Fermi velocity. 

The correlation length critical exponent is then 

"'^'^ . + 0(e"). (11) 



2 {7-S){9 + S) 



and 



Vb 



7-S 



(12) 



with r]i, as the usual OP's anomalous dimension. 

The fermion propagator at the critical point behaves as 



dimension is 



' + k'^)^^ where the fermionic anomalous 



77/ 



2(7-5') 



(13) 



The scaling then implies^ that the residue of the Dirac 
quasiparticle's pole vanishes near the critical point as a 
power-law. 



(14) 



There are therefore no sharp fermionic excitations right 
at the critical point. 

Both the bosonic and the fermionic anomalous dimen- 
sions, being proportional to the critical point values of 
the Yukawa coupling, are finite already to the leading 
order in e. This implies that the exponents have rather 
different values from the usual mean-field. Since there 
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does not seem to be a dangerously irrelevant coupling in 
the problem the hyperscaling should hold and we may 
obtain the remaining exponents from the usual scaling 
relationsj^i This way we find 



4(3 + 5') ^, 2^ 



(15) 



<5=4i^ = 3-i±4. + 0(a. (16) 



d-1 



Vb 



7-5 



For a large number of Dirac fields N one similarly finds ;^ 
ly = 1 + 0{1/N) and rjb = 1 + 0{1/N), and thus 7 = 
1 + 0(1 /N), and (5 = 2 + 0{1/N). Various calculations 
on the Gross-Neveu models in the pasl3 also suggested 
that the infinite-N results are often good estimates of the 
exponents' actual values. 



IV. DISCUSSION 



from the insulating side both the masses of the Higgs 
OP and of the Dirac fermions approach zero, with the 
universal ratio 



TJif gl 9 



16 



5 



(17) 



analogous to the universal Ginzburg-Landau parameter 
at the superconducting critical point. The emerging 
Lorentz invariance in the critical region implies that a) 
the dynamical critical exponent z = 1, so that the Fermi 
velocity is not critical)^ b) the (non-universal) velocity of 
bosonic excitations in the critical region approaches the 
(non-universal) Fermi velocity. For the triplet OP this 
means that the spin- waves on the insulating side will have 
the same velocity as the Dirac fermions on the metallic 
side, when both are near criticality. As a consequence, 
the ratio of the specific heats on the metallic and insu- 
lating side in the critical region approaches a universal 
value 



The critical exponents have recently been computed 
numerically for two species of Dirac fermions interact- 
ing via ^ 1/r interactions. Whereas only the stan- 
dard power-laws and no essential singularit y^ ^1^^ were 
observed)^ the critical exponents 7 = 1, and 8 ~ 2.3, 
particularly when expressed in terms oi v = 0.85 and 
rib = 0.82 appear distinctly not to be of the usual mean- 
field variety. They do seem close to our one-loop re- 
sults, however, particularly for the triplet OP, where for 
e = 1 we find v « 0.88 and rjb k, 0.80. The correlation- 
length exponent is also very close to the previous two- 
loop result^Si. for the singlet OP. These are only crude 
estimates, but we believe that the observation of the val- 
ues of V and r]b near unity should be taken as a sign of 
the Gross-Neveu- Yukawa fermionic criticality. Although 
the strength of 1/r interaction may be used as the tun- 
ing parameter, both the present e-, earlie r^i^^i^^ 1/A^- 
expansions, and the results near 1-1-1 dimensions'^ sug- 
gest this to be ultimately an irrelevant coupling. The 
Mott transition in graphene may be analogous to the sit- 
uation in the QEDi^ where the CS breaking transition 
brought by the increase of the electromagnetic charge 
appears to be described by the Nambu-Jona-Lasinio La- 
grangian with only short-range interactions. This is one 
way to phrase the triviality of the continuum limit of 
QED4^ 

Let us list some further results specific to the Mott crit- 
icality in graphene. As the critical point is approached 



The numerator (denominator) corresponds to the spe- 
cific heat of the eight (two) free relativistic fermions 
(Goldstone bosons) in d dimensions. Assuming that the 
Lorentz invariance is still emergent at the critical point 
in d = 3^1 the same ratio is exactly 3 at the transition 
into the spin-density-wave phase in the Hubbard model 
on honeycomb lattice. 

In conclusion, we presented an effective theory of the 
Mott transition on honeycomb lattice that becomes solv- 
able near three spatial dimensions. We argued that the 
quantum critical point in graphene should be relativisti- 
cally invariant, and possibly logarithmically trivial. An 
interesting consequence of our theory is the existence of 
a simple universal ratio of the specific heats of the two 
phases near the transition. This particular prediction 
should hopefully become testable in the future numerical 
and experimental studies. 
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